Divergence of the effective mass near a density wave instability in a MOSFET system. 
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We study the renormalization of the Fermi-liquid parameters in the vicinity of a density wave 
quantum phase transition, which should occur in MOSFET systems at low densities. First, using a 
perturbative RPA treatment of fluctuations, we calculate the electronic self-energy and show that 
the effective mass diverges at the density wave transition point. Second, we go beyond perturbation 
theory, making use of the exact Pitaevskii identities. Within this exact analysis, we also find 
a divergence of the effective mass, which occurs at higher densities in the fluctuation region, as 
compared to the perturbation theory. This result signals the break-down of conventional Fermi- 
liquid description in the vicinity of the transition point. The divergence of the effective mass gives 
rise to a singular behavior of the electronic compressibility. We suggest that the experimentally 
observed enhancement of the effective mass is a precursor to a second order thermodynamic phase 
transition into a glassy density wave state. 
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I. INTRODUCTION 

Novel experiments in two-dimensional high-mobility 
electron systems have provided a great variety of re- 
markable results which do not fit into the framework of 
the conventional theory of fermion systems i The most 
prominent effect is indeed the metal-insulator crossover 
observed in such systems, which clearly contradicts to 
the scaling theory of localization 2 and a lot of attempts 
have been made to reveal the underlying nature of 
the phenomenon. Apart from the unexpected metal- 
lic regime, there is a number of other surprising ex- 
perimental results, such as the sharply increasing effec- 
tive mass at low electron densities^*^ in silicon metal- 
oxide-semiconductor field-effect transistors (MOSFETs). 
There are also experiments clearly showing some unusual 
behavior in thermodynamic properties such as the local 
electronic compressibility, 6 suggesting the existence of a 
second order thermodynamic phase transition in the sys- 
tem, which, if indeed exists, should be related to the 
quantum crossover observed in the transport behavior. 

A successful theory describing the peculiar physics of 
the high-mobility electron systems must be based on the 
correct choice of the quantum ground state of the sys- 
tem. It is clear now that traditional paramagnetic Lan- 
dau Fermi liquid theory in its perturbative form does not 
provide an adequate explanation of the observed physics. 
One is forced to assume the existence of a density-driven 
quantum phase transition(s) in the two-dimensional elec- 
tron system. 

The competition between the potential and kinetic en- 
ergies of a pure electron system yields the first order 
phase transition into the Wigner solid phase at very low 
densities. This fact was realized a very long time agOj£ 
but the theoretical description of the liquid-solid phase 
transition is still missing and only numerical studies are 
available. These numerical studies derive extremely high 
values of the Wigner transition RPA parameter: r s ss 110 



and r s ~ 37 in thread and two^ dimensions correspond- 
ingly. There are also numerical works^Siii*^ which sug- 
gest that in between Wigner crystal and usual Lan- 
dau Fermi-gas there may be an intermediate phase, the 
charge— or spin-density wavei^ phase, in which a sta- 
ble charge (spin) density modulation exists in the ground 
state with some wave- vector Q < 2pp: 



M ! 



(r) = n [1 + A cos (Qr + %)] 



(1) 



where for the charge-density wave: <5j = <5j and for the 
spin-density wave: 6^— 5± =tt. Parameter A serves as the 
order parameter for the phase transition. Usually, a den- 
sity wave transition is associated with the Peierls insta- 
bility in effectively one-dimensional systems^2ii4 which 
occurs due to the strong electron-phonon coupling. Nev- 
ertheless, such a state may appear in higher dimen- 
sional electron systems within a jellium-like model, as 
was first pointed out in Ref. [l^. At low electron den- 
sities, there is a region where a density wave state is 
energetically more favorable compared both to the uni- 
form density distribution corresponding to usual Landau 
Fermi gas and to Wigner crystal i&iLH Numerical simu- 
lations of Swierkowski et alw— showed that in double- 
layer and multi-layer systems, the transition into the 
charge-density wave state occurs at r s ~ 10 — 20, i.e. 
much earlier than the Wigner crystal phase transition is 
expected. The charge-density wave state in an isotropic 
system is possible only if the neutralizing "jellium" is 
polarizable, while the spin-density wave should occur for 
any elasticity of the positive background. In a MOSFET 
system, the only possible candidate for such a polariz- 
able background is the charge distributed on the metal- 
lic gates, which is separated from the two-dimensional 
system of interest by an oxide dielectric barrier, usually 
Si02, of some finite thickness d. Indeed, if the thick- 
ness is very large compared to the typical inter-electron 
distance, the "jellium" on the gates does not "see" any 
changes of the electron density at the length-scale of 1 /p-p 
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and any charge-density modulations should become un- 
stable. However, if the oxide layer is thin enough, there 
should exist a critical density at which a phase transition 
into the charge-density phase occurs. This corresponds to 
the experiments of Tsui et al.^ where a divergence of the 
effective mass was reported for a thin-oxide MOSFET. In 
the majority of other experiments, the thickness of the 
oxide is about d ~ 10 3 A, which is much greater than the 
typical inter-electron distance, and the transition into the 
spin-density wave phase is anticipated. Thus, indepen- 
dently of the geometry of a two-dimensional structure, 
a density-wave phase may exist separating the paramag- 
netic Fermi liquid and Wigner solid phases. 

The subject matter of the present paper is to explore 
the behavior of the Landau Fermi-liquid theory param- 
eters as one approaches a density wave transition. Our 
main finding is that the Fermi-liquid theory apparently 
breaks down in vicinity of the transition point as the ef- 
fective mass diverges, which leads to the instability of the 
usual Fermi-step density distribution even at T — 0. We 
suggest that the observed enhancement of the effective 
mass^^ is a precursor of a density wave transition in 
the two-dimensional systems. 

Our paper is structured as follows: In Sec. II, we study 
the renormalization of the quasiparticle Z-factor and the 
effective mass within the self-consistent RPA approxima- 
tion in the vicinity of the charge-density wave transition. 
This method of treating fluctuations is in the spirit of the 
classical paper of Doniach and EnglebergjiSi who consid- 
ered the ferromagnetic transition case, and also similar 
to Ref. ^3 j where the self-energy function was studied 
near the antiferromagnetic transition in connection with 
the physics of high-T c cuprates. First, we derive gen- 
eral expressions for the self-energy function. Second, we 
discuss the structure of the appropriate RPA density- 
wave fluctuation propagator. The actual calculations are 
done for the charge density wave case and it is shown 
that within the RPA approach the quasiparticle Z-factor 
vanishes at criticality (c. /., Ref. and the effective 
mass diverges as 1/Z at the critical point. In Sec. Ill, 
we study the fluctuation effects going beyond conven- 
tional RPA perturbation scheme, by using some exact 
formulae from the Landau Fermi-liquid theory. Namely, 
on the basis of the Pitaevskii identities, we prove that 
for a density-wave phase transition (and more generally, 
for any phase transition, which breaks the translational 
invariance of an otherwise isotropic system) the renor- 
malization of the derivative {dYi/duf) is never singular 
and the Z-factor formally never vanishes as distinct of 
the RPA result. The effective mass does diverge in the 
fluctuation region near the phase transition, but not at 
the criticality itself, again, opposed to the RPA conclu- 
sion. At this point, the group velocity of quasiparticles 
changes sign and the conventional Fermi-liquid theory 
breaks downiiSiiS Moreover, at this point some thermo- 
dynamic quantities, such as the electronic compressibil- 
ity, experience a singular behavior. In Sec. IV, which 
is the conclusive part, we summarize our findings and 



briefly discuss our understanding of possible physics in 
the density- wave phase. 



II. RPA-TREATMENT OF THE 
DENSITY- WAVE FLUCTUATIONS 

A. Fermi-liquid theory distilled 

The effective Fermi-liquid parameters follow from the 
quasiparticle Green function, defined as usuali^S 



G c (*i,r i; t 2 ,r 2 ) = -i (Tip (ti,^)^ (*a,r 2 )) 



(2) 



where T is the time-ordering operator and i/j and rjv are 
electron field operators. In what follows, we will use sub- 
script "c" to denote a time-ordered quantity and sub- 
scripts "R" and "A" to denote a retarded and advanced 
quantities correspondingly. The Green function @ is in 
turn determined by the self-energy function S (e, p) as 
follows: 



Gc (e,p) 



1 



e-e o) (P)-S c (e,p) 



(3) 



In the vicinity of the Fermi-momentum, Eq. (|3J) can be 
re-written as 



G c (e,p) 



£-£(p)±«7 (e,p)' 



(4) 



where Z € [0 , 1] is the jump in the quasiparticle distri- 
bution function at the Fermi level, which is determined 
by the following identity: 



d 

— ReS R (s,pf) 



(5) 



e=0 



and £ (p) = (|p| — pf)j>f/w* is the quasiparticle spec- 
trum, with the effective mass defined via the following 
relation 



m 



d 



Re£ R (0,p) 



(6) 



The parameter 7 = ZImE c (e,pf) determines the quasi- 
particle life-time and should be small compared to the 
energy of a quasiparticle, for the latter to be well-defined. 
Under this assumption, we can write the retarded Green 
function in the following form, which shall be used in 
further actual calculations: 



Gr (e,p) = P. 



Z 



£(p) 



- mZ5 [e - £ (p)] , (7) 



where symbol P stands for the principal value of the 
corresponding integral. 

Within the perturbative treatment of fluctuations, 
only the first order diagram (containing one fluctua- 
tion propagator) for the self-energy is retained. The 
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corresponding analytically continued retarded self-energy 
function reads (for e > 0)t2i 



(e,p) 



1 



+ 00 

2 " ' duj 



d 2 q 



x < Im Gr (oj + e, p — q) Dr. (— cj, q) tanh 



oj + e 



2T 



+ G R (oj + e, p - q) Im D R (uj, q) coth 



2T 



(8) 



The interaction D(u>,q) changes the properties of the 
system and renormalizes the Fermi-liquid theory param- 
eters. In the case under discussion, the density wave 
fluctuation propagator will be used as the appropriate 
effective interaction term. 

From now on, we will consider zero temperature case 
only, so that the tangent and cotangent factors in Eq. © 
reduce to theta-functions. Moreover, using Eqs. J7J and 
iJSJ , we can evaluate the integral over the directions of q 
exactly (provided that the interaction is isotropic) and 
obtain the following expression for the real part of the 
self-energy, which is quite involved: 



ReE R (e,p) = E 1 (p,e) + E a (p,e), 



(9) 



where 



Si (e,p) 



Zv 

7T 



and 

£2 (e,p) 



qi («) 

dui sgn (oj + e) J dq 
q Re Dr (oj, q) 

x Vb 2 - <ii M [qi W) 1 2 } 

00 f gi(Sl) 00 

— j doj sgn oj\ J dqq + J 

-00 [ g 3 (fi) 

Im D K (oj,q) 



(10) 



dq 



,(n) 



where 

qi,2 (A) = V2pf 



2El 



T 1 



1/2 



1/2 



(12) 

and f7 = (oj + £ — £ p ), with ^ being the density of states 
at the Fermi-line and Ep the Fermi energy. 

Using Eqs. JJjJ, (|10|) . and one can prove the fol- 

lowing useful identity: 

_ 2Zi/ ^ dqD(0,q) 

P = Pv = ~ 
e = 



_9 _a_ 

9e 



Er (e,p) 



V( 2 Pf) 2 - 9 2 



(13) 

Let us note, that the right-hand side of Eq. Ijl^jl is very 
large, provided that the effective propagator is singular 
in the static limit. 



B. Charge-density wave fluctuation propagator. 

Now let us consider a thin-oxide MOSFET in which 
the two-dimensional electron system is separated from 
the metallic gates by a barrier of width d, which we sup- 
pose to be not too large compared to the typical inter- 
electron distance. In this case, the charge on the metallic 
gates serves as an elastic (not rigid) neutralizing back- 
ground, which can compensate the energy cost produced 
by inhomogeneities in the electron density distribution 
in the two-dimensional electron system. At sufficiently 
low electron densities, one anticipates a transition into 
the charge-density wave state Q which should be ener- 
getically more favorable than the homogeneous distribu- 
tion. Below the critical density n c the initial Galilean 
invariance of the system is broken. However, above the 
transition the system is still isotropic and only fluctuat- 
ing density waves can exist with no preferable direction. 
The appropriate static fluctuation propagator, quite gen- 
erally has the form: 

D(uj = 0,q)^-— —2, (14) 



(q-QY 



where £ is the charge-density wave coherence length, 
which diverges at the transition and Q ~ 2pp. The static 
part of the propagator l|14|) has a typical Lorentzian form 
(see, e g., Ref. |22|) and, as we noted, should be still 
isotropic above the transition as the Galilean invariance 
is preserved until the transition point: n = n c . The 
transition itself is thermodynamic in nature and is all 
about the sign of the difference between the two ener- 
gies of the homogeneous distribution and distribution Q . 
The RPA philosophy assumes that the dynamic part of 
the fluctuation propagator, which is not crucial for the 
transition itself, can be obtained by summing up a series 
of diagrams containing the RPA polarizability bubbles, 
which generate the required frequency dependence. Tak- 
ing into account these considerations, we write the RPA 
charge-density wave fluctuation propagator in the follow- 
ing form: 



A 



D(oj,q) = 



vo a + [1 - q/QY + iP (oj, q) ' 



(15) 



where a ~ (n — n c ) /n c is the Ginzburg-Landau-like co- 
efficient, which changes sign at the transition point, vq 
is the bare density of states at the Fermi-line and the 
imaginary part iP has the form 

P(oj,q) = -ImU(oj,q) 

where II (oj, q) is the polarizability bubble and A and A 
are some dimensionless constants, which may depend on 
r s only. The polarizability is defined as usual (here, we 
use Matsubara notations): 

/^ 2 k 
- — -2<? (e„, k) G (e n + oj m , k + q) . 
(27T) 



(16) 
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For a two-dimensional electron gas at T = 0, this quan- 
tity was calculated by SternS^. We do not need its exact 
expression. Let us only emphasize, that the case Q ~ 2p F 
may be different from the case Q = 2p F , since the polar- 
izability is a non-analytic function of frequency at q = 
2p F and the mixture of the two non-analyticities makes 
the problem technically more cumbersome, (for Stern's 
polarizability: 23 Imll(a;, 2pp) ~ vM s S ncJ )- Quite gen- 
erally, Q < ppr"^ hi which case we can expand the 
propagator at low frequencies to obtain: 



vZ* ( — , as uj -► 0, 



where ujq is a constant, which does not depend on the 
closeness to the critical point, Z is the quasiparticle Z- 
factor, and v — m*/ (2tt) is the density of states, both to 
be found. 

In the case of a spin-density wave transition, which 
may occur in a MOSFET with a large oxide barrier, the 
propagator should have the form similar to 1)14(1 but with 
a non-trivial spin structure: 



Z- 2 + (q-Qy+tP(uJ,q)' 



(17) 



In what follows, we will do explicit calculations for the 
case of the charge-density wave only. Generalization to 
the spin-density wave case is straightforward and the 
qualitative results are identical. 



C. Renormalized Z-factor and effective mass 

Now, we are at the position to calculate the renor- 
malization of the Fermi liquid parameters due to charge- 
density fluctuations within the RPA-approach. Identity 
((13)1 allows us to calculate only one derivative, the other 
will be automatically extracted from Eq. 1)13)1 . Let us 
focus on the e-derivative, which determines the renor- 
malization of the quasiparticle Z-factor. From Eq. I)l(jp. 
we have 



—Ei £,p = — 

as 7r 



91 M 



dw sgn uj 



dqq 



d 

x — Re D R (w 

oe 



e,q), (18) 



where functions qi.2{^) are defined in Eq. (JT2J. Let 
us consider the region in the immediate vicinity of the 
transition only so that a <C (2p F — Q) /2pp. Then, the 
integrals in Eq. 1)10)1 and 1)11)1 are determined by low fre- 
quencies u) ~ a/Z 2 — > 0. Using the following properties 
of any physical propagator 



and 



Im £> R (w, q) = - Im D R (-u>, q) 



Re D R (ui,q) = Re D R (-u), q) . 



one can see that expression 1)10)1 is identical to the right- 
hand side of Eq. 1(13(1 . Eq. ((11(1 is exactly zero as the prin- 
cipal value of the corresponding integral. In the leading 
order, we obtain the following results: 



d_ 

de 



ReE (s,pf) 



2A 



e=0 



and 



— ReE(0,p) 



i - (Q/2 PF y 



= o (a) —>■ 0, 



(19) 



(20) 



P=PF 



where we remind that a = [n — n c ) /n c is a small devia- 
tion from the charge-density wave transition point. 

Let us mention that besides the singular fluctuation 
effects we studied so far there are other contributions 
coming from the q — > channel. Such corrections have 
been considered previously by many authors (see e.g., 
Ref. |24|) and it was found that they do increase the ef- 
fective mass by a factor of two or so (Ting et al. 24 predict 
m*/m sa 2.0 for r s = 5) but do not lead to any divergent 
behavior. Let us denote the corresponding renormalized 
effective mass as rh and the quasiparticle Z-factor as Z. 
In these notations we obtain the following self-consistent 
equation for the quasiparticle Z-factor within the RPA 
approach [we consider here the limit a -C (2pp/Q — 1) ]: 



2A 



Z 2 (a) Z(a) 



1 - (Q/2p F ) 



Z 



1 = 0. 



(21) 



This quadratic equation can be easily solved exactly, but 
we are mostly interested in the behavior in the closest 
vicinity of the transition, which leads to (c. /., Ref. [l7|): 



Z(a -> 0) = 



Q_ 

2 PF 



a 1 / 4 oc (n — n c ) 



1/4 



(22) 



Thus, the effective mass diverges as 1/Z, and we have 



1 



' («) (n - n c ) 



jn, Q<2p F . 



(23) 



We see that within the RPA-treatment of charge-density 
wave fluctuations, the effective mass diverges at critical- 
ity, while the Z-factor vanishes. 

Let us mention that the case Q — 2p F is somewhat 
pathological, since the pole of the propagator coincides 
with the point of the non-analyticity of the polarizabil- 
ity function and also with the point where the inverse 
square-root term appearing in the two-dimensional inte- 
grals lfTU|l and (fTT|) diverges. To obtain the correct nu- 
merical factors and to avoid unphysical divergences, one 
should use the exact form of the propagator in the hole 
range of the (w, q)-space. If we are not interested in ex- 
act numbers, we can easily estimate the scaling law of the 
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Z-factor near the transition for the Q = 2p F case within 
the RPA approximation: 



Z Q=2PF (a) ~ A 1 / 2 ^ 8 
Thus, the effective mass diverges as 

f 



(n - n c ) 



3/8 ■ 



Q = 2p F 



(24) 



(25) 



Let us emphasize that results <|22l — - I25|) take into account 
only one diagram, which alone yields very divergent re- 
sults. 



III. NON-PERTURBATIVE ANALYSIS. 

Although, the RPA treatment is widely-accepted and 
often leads to reasonable results, we believe that in some 
cases there are serious doubts concerning its reliability. 
Unlike in the clcctron-phonon problem, where the Migdal 
theorem holds and circumvents the necessity to study 
complex vertex corrections, some other problems (such 
as, e. g.: paramagnon coupling near a magnetic phase 
transition^&ii charge density wave transition, spin den- 
sity wave transition, etc.) require correct account for the 
vertex corrections, which are not small in the vicinity 
of the critical point and the Migdal theorem is appar- 
ently violated. The exact account for the vertex correc- 
tions seems an insurmountable problem at the moment 
and the usual practice is to cover the underlying diffi- 
culties by saying that the domain of applicability of the 
RPA-like treatment may be extended and that the per- 
turbative approach should qualitatively explain the key 
physics of the transition. In the present section, we prove 
that such an arguing is quite dangerous and may lead 
to some qualitatively incorrect conclusions. Fortunately, 
the density-wave transition case may uncover the under- 
lying problem, as one can obtain some very important re- 
sults based on very general grounds. As we show below, 
in the vicinity of any phase-transition, which breaks the 
Galilean invariance of an initially isotropic system, some 
important effects are hidden in higher order diagrams. 
The main qualitative result we are going to derive in the 
present section is that the effective mass diverges before 
the phase transition occurs and this divergence leads to 
the break-down of the Fermi-liquid theory. Apart from 
the perturbative result, predicting Z = at the critical- 
ity, we show that the Z-factor formally always remains 
finite within the non-perturbative treatment. 

Paradoxally, we will use the essence of the Fermi-liquid 
theory construction to justify its break-down. Let us in- 
troduce such a standard element as the irreducible four- 
vertex function^ T(j>i;p2\pi + k\p 2 + fc), where we in- 
troduce notations pi = (si,Pi) for brevity, with pi and 
P2 being the momenta of incoming and outcoming par- 
ticles correspondingly. In the limit k = (w,k) — > the 
four-vertex has a singular structure and we, following 



Ref. [2l| introduce the following standard function: 

r w (pi;p 2 )= hm T{p 1 ;p 2 \ Pl +k;p 2 + k). (26) 

|k|/w->0 

We also decompose the product of two Green functions 
involved in the calculation of the four-vertex T into a 
singular and regular parts i^i 



G(p)G( P + k) 



2niZ 2 vk 
v lu — vk 



8(e)8 (|p| -pf) + 5rog(p), 
(27) 

where <?rcg(p) is a regular function in the limit k, lu — > 0. 
One of the central identities of the formal Landau Fermi 
liquid theory construction reads^i 



dG' 1 (p) 
de 



1 

~Z 



1 + - T»(p,p') 9les (p')—P 
£ J i (2tt) 



(28) 

The other identities, which we refer to as the Pitaevskii 
identities^ can be obtained from Eq. (|28|) by consider- 
ing symmetries of the initial Hamiltonian (in our case 
isotropic electron system with Coulomb interactions be- 
tween electrons). The Galilean invariance, which is bro- 
ken in a density wave state, but is preserved above the 
phase transition yields the second Pitaevskii identity 
which readstSii2^ 



1 



d 3 p' 



p 9 -^ 1 = p + i J r- tort pWpO t^s (29) 

From Eqs. if^Sl) and the following relations follow 

immediately: 



<9E (s,pf) 



0e 



r w (p;p')g r e B (p') 



d 3 P ' 

i(2nf 



(30) 



and 



^ (e, PF ) 



0s 



, t\ (ppO 



Pi 



ip') 



d 3 P ' 

i (2tt) : 



(31) 

In the case of a phase transition, preserving the Galilean 
invariance of the system, the four- vertex function should 
be peaked at |p — p'| = and the relations l|30|l and l|31() 
become identical in the leading order, since (pp') = pp. 
In the case of a second order phase transition, which does 
break the translational symmetry (of which a density- 
wave phase transition is an excellent example indeed) 
we anticipate the four-vertex to diverge at the transition 
point at some finite wave- vector Q. I. e., in the fluctu- 
ation region, the major contribution should comes from 
the domain |p — p'| = Q. Therefore, relations l|30|l and 
(|31|l lead to an obvious contradiction, which is especially 
pronounced in the case Q = 2pp, when (pp') jp\ = — 1. 
The only possible resolution of the paradox is to accept 
that 



9£ (e,p F ) 



0e 







(32) 



e=0 
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in the leading order in the closeness to the phase transi- 
tion. We see that this exact result is in some sense op- 
posite to the one obtained within the perturbative treat- 
ment [see Eq. <|22|) ]. which proves that the RPA-approach 
contains a serious deficiency. Obviously, the singularity 
can not just disappear from the theory, since the general 
Landau relation-^ ensures that at least one of the deriva- 
tives ((E/de) and {cE/dp) is singular [c. /., Eq. lfT3|) ]. 
The relation reads 



d_ _d_ 

de d£ p 



S(e,p) 



P = PF 

e = 



d 2 p 

(27rp F )' 



r/(p,p') 



dnlp') 
X ^ P ' 



(33) 



where we introduced the full scattering amplitude (Lan- 
dau function), which is the static four- vertex function: 

/(p,p')=r" (0,p;0,p'). 

In the perturbative approach, this singularity arose in the 
e-derivative of the self-energy. We see, that within the 
very general non-perturbative analysis the singularity is 
"transferred" into the p-dcrivative, which automatically 
means that the quasiparticle Z-factor never vanishes and 
the only quantity which acquires singular contributions 
and eventually diverges is the effective mass. The equa- 
tion determining effective mass renormalization can be 
derived from Eq. l(2"%|l and the gauge invariance of the 
system, and in the most convenient form can be written 
as: 



dp 



- f Z(p-p') 



dn(p') d 2 p' 
dp' (2tt) 2 



(34) 



where n(p) is the electron distribution, which does not 
have to be the Fermi-distribution function, in general. 
£ p [n] is the electron spectrum, which in general is a com- 
plex functional of the distribution functionii£*i2*Si Ide- 
ally, the electron distribution function should be self- 
consistently determined from the Landau energy func- 
tional, as its absolute minimum and must be subject to 
the stability constraint d£, p [n]/dp > (i. e., to* > 0). 
However, as long as the group velocity of the quasipar- 
ticles (the inverse effective mass) is positive, the familiar 
Fermi-step solution is the correct distribution function at 
T = 0. 

The scattering amplitude / is a singular quantity in 
the vicinity of the transition and quite generally can be 
written as: 



/(P-P') = --|« + 



Ip-p'I-Q 



(35) 



where again Vq = m/ (27r) is the "bare" density of states 
at the Fermi line, a oc (n — n c ) is a small parameter, 
which changes sign at the transition, A is a dimensionless 
constant, which may depend on r s only, and (2tt/Q) is the 



spatial period of the density wave modulations. Eq. I|34() 
leads to the following formula for the effective mass: 



1 1 2A 

TO TO* TO7T 



sin 1 

Q 2 



cosx^X- 

After elementary integration we obtain the following re- 
sult for \/2pp < Q < 2pp: 



1 - 2A 



\Q 2 /(2 P - 



1 1 



tJ(2 Pf /Q) 2 -1 



(37) 



and for Q = 2p-p'. 



to* = m 



1 - 2Aa- 3 / 4 



(38) 



We see that in the vicinity of the transition a — > 0, the 
effective mass diverges, independently on the value of A, 
and formally changes sign after this point in the fluctua- 
tion region. This immediately signals that the Fermi-step 
distribution function is not stabl o 18 ! 19 and conventional 
Fermi liquid theory breaks down. 

Another quantity, which can be easily estimated from 
the general considerations of the Landau Fermi liquid 
theory is the electronic compressibility, which can be ex- 
pressed as2& 



\Sn) to* ( 



1 + */(*)) 



where v — m*/ (2n) is the density of states and / (%) is 
a singular quantity in the vicinity of the phase transition 
[see Eq. I)35[l]. Therefore, we conclude that the electronic 
compressibility \dji/dn\ must show a divergent behavior 
in the vicinity of the phase transition (as to* — * oo). 
Clearly, other thermodynamic quantities such as specific 
heat, susceptibility (g- factor), etc. may also acquire a 
singular behavior. A detailed analysis of thermodynam- 
ics in the vicinity of the charge- and spin-density wave 
transitions will be reported elsewhere. 



IV. CONCLUSION 

Summarizing, we suggest that the metal-insulator 
transition may be actually a signature of a thermody- 
namic phase transition in a dilute two-dimensional elec- 
tron system. Namely, the charge- or spin-density wave 
transition. Which of the two phases prevails is a tenta- 
tive question and the answer strongly depends on the 
geometry of the system and can be determined only 
on the basis of an accurate numerical analysis compar- 
ing the energies of the two states. Experimentally, the 
spin-density wave case may be distinguished from the 
charge-density wave one by measuring the electron g- 
factor which should be singular in the former scenario and 
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unremarkable in the latter. Our expectation is that in 
thin-oxide MOSFET's^ the charge-density wave should 
be more energetically favorable, while the spin-density 
wave should take over if the metallic gates are located 
far from the two-dimensional electron gas. In both sce- 
narios, we proved that the density-wave fluctuations lead 
to the divergence of the electron effective mass and to 
a singular behavior of the electron compressibility. Both 
effects have been observed experimentally^i^ The den- 
sity corresponding to m» = oo is the point, where the 
Fermi-liquid description in its conventional form appar- 
ently breaks down. In the immediate vicinity of the den- 
sity wave transition, the problem directly maps on the 
model of so-called fermion condensation introduced ear- 
lier by Khodel and ShaginyanA 8 " and considered later by 
NozieresiS and also by Volovik^ The issue of such a non- 
Fermi-liquid state is far from its closure, but the main 
consequence is that the Fermi-distribution function can 
not be a stable solution of the Landau energy functional 
if it leads to a negative group velocity of the quasipar- 
ticles. Possible model distributions were suggested; 18 i 19 
all characterized by a very high density of states at the 
chemical potential v{p) (with the singularity smeared out 
only by temperature) , which may correspond to the sharp 
increase of the conductivity. 

Below the phase transition, i. e., in a density wave 
state, disorder should play the key role leading to the 
destruction of the long-range order and to the pinning 
of density waves (according to the general theorem due 
to A. I. Larkin 28 ). The glassy nature of the density 
wave state should reveal itself via a noisy behavior of 
linear response quantities (such as conductivity), which 
may show some real-time glassy behavior. The physics 
here should be quite similar to the physics of supercon- 
ducting vortices in the presence of disorder piii^SSi Let 



us mention experiments of Jaroszyhski et al£SL and of 
Bogdanovich and Popovicj^i in which the low-frequency 
resistance noise was measured clearly showing the glassy 
freezing of the electronic system. 

Finally, let us emphasize an important technical dis- 
covery concerning the reliability of a perturbative RPA- 
like treatment of strong fluctuations and the impor- 
tance of vertex corrections. The problem is quite gen- 
eral and is related, in particular, to the physics of fluc- 
tuation phenomena near the antiferromagnetic phase 
transition in high-T c cupratestf We have shown that 
the vertex correction, which are usually excluded from 
consideration^&i! yield qualitatively important changes 
and ensure that the quasiparticle Z-factor does not van- 
ish. Our current understanding, based on studying a 
model example allowing parquet technique treatment (to 
be published elsewhere), is that the effect of the vertex 
corrections is mostly the renormalization (suppression) 
of the imaginary part of the initial RPA-like propagator. 
This result seems quite reasonable, since the imaginary 
part of the propagator is related to the life-time of fluc- 
tuations, which should become infinite at the transition 
point. 
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